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Abstract 



We report the resuhs of a study of two-body B meson decays to the complete 
set of Kit, vrvr, and KK final states. The study is performed on a data sample 
of 31.7 lb 0.3 million BB events recorded on the T(4S) resonance by the Belle 
experiment at KEKB. We observe significant signals in all Ktt final states and 
in the tt'^'tt" and 7r"'"7r'' final states. We set limits on the vr'^vr" and KK final 
states. A search is performed for partial-rate asymmetries between conjugate 
states for fiavor-specific final states. 
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I. INTRODUCTION 



There exists a gold mine of weak and hadronic physics in two-body B meson decays to 
the Kn, vrvr, and KK final states. Indeed, if one assumes unitarity of the quark mixing 
matrix [|l|, these modes contain enough information to measure all angles of the Unitarity 
Triangle f^. Methods to extract weak-sector physics from these decays are complicated by 
hadronic uncertainties. However, if enough final states are measured, we will have sufficient 
information to constrain the sizes of hadronic amplitudes and strong phases, a necessity in 
disentangling the unitarity angles from measurements of fiavor asymmetries and the relative 
size of the partial widths among these modes PHIO . 



We have previously reported measurements of, or limits on, the branching fractions of B 
mesons to the Kn, vrvr, and KK final states excluding the tt^tt" and K^K^ final states |TT 



as well as a search for charge asymmetries in the fiavor-specific Kn final states |jT2[, based 
on a data sample of 11.1 million BB events. The results presented here include the previous 
data and supersede all previous results. Similar studies have been performed by other 
experiments [p!^-pO| • 



Here, we present measurements of, or limits on, the branching fractions of B mesons to 
the Kn, vrvr, and KK final states including all combinations of charged and neutral kaons 
and pions. We refer to these final states collectively as -B ^ hh, including charge conjugate 
states unless explicitly stated. For final states where the charge of the kaon or pion specifies 
the fiavor of the parent B meson, known as fiavor-specific final states, we present limits on 
the partial-rate asymmetries defined as 

N{B ^ f) - N{B ^ f) 

where B represents either a B^ or B~^ meson, / represents a fiavor-specific final state, and 
B and / are their conjugates. 



II. APPARATUS AND DATA SET 



The analysis is based on data taken by the Belle detector [|2T]] at the KEKB e^e storage 
ring p2|. The data set consists of 29.1 fb~^ on the T(4S) resonance corresponding to 



31.7 dz 0.3 million BB events. An off-resonance data set of 4.4 fb~^ was taken 60 MeV below 
the T(4S) resonance to perform systematic studies of the continuum e^e^ qq background 
where q is either a u, d, s, or c quark. KEKB collides 8 GeV electrons and 3.5 GeV positrons 
that are stored in separate rings, producing an T(4S) system that is boosted by 7/3 = 0.425 
along the beam axis. In this analysis, all variables are calculated in the center-of-mass frame 
of the electron and positron beams unless explicitly stated. 
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The Belle detector is a general purpose magnetic spectrometer with a 1.5 T axial magnetic 
field. Charged tracks are reconstructed using a 50 layer central drift chamber (CDC) and 
a 3 layer double-sided Silicon vertex vetector (SVD). Candidate electrons and photons are 
identified using an 8736 crystal CsI(T/) calorimeter (ECL) inside the magnet. Muon and K^^ 
candidates are identified using resistive plate chambers embedded in the iron magnetic fiux 
return (KLM). Hadron and auxiliary lepton identification is provided by an array of 1188 
Silica aerogel Cerenkov threshold counters (ACC) and a barrel of 128 time-of-fiight (TOF) 
plastic scintillator modules. 



III. EVENT RECONSTRUCTION 

Event triggers based on fast signals from the CDC, ECL, TOF, and KLM ETI]. Hadronic 



events are selected using event multiplicity and total energy variables p3[- For signal events 
that pass all B ^ hh selection criteria, the triggering and hadronic event selection efficiencies 
range from 99% for — > h^h~ modes to 76% for the B^ — »• 7r°7r° final state. 

Charged vr and K mesons are identified by their energy loss {dE/dx) in the CDC and 
their Cerenkov light yield in the ACC. For each hypothesis {K or vr), the dE/dx and ACC 
probability density functions are combined to form likelihoods, Ck and C^^. K and vr mesons 
are distinguished by a cut on the likelihood ratio Ck/{Ck + Ct,). A similar likelihood 
ratio including calorimeter information is used to identify electrons. All charged tracks 
that originate from the interaction point and are not positively identified as electrons are 
considered as kaon or pion candidates. 

Candidate mesons are reconstructed using pairs of oppositely charged tracks that 
have an invariant mass in the range 480 MeV< m(7r+7r^) < 516 MeV p4|. The candidate 
must have a displaced vertex and fiight direction consistent with a originating from the 
interaction point. Candidate 7r° mesons are formed from pairs of photons with an invariant 
mass in the range 114 MeV< "1(77) < 150 MeV and -E'y(iab) > 70 MeV. The tt+tt^ and 
77 mass spectra are shown in Fig. [l| for B^ K^h^ and /i+7r° candidates in the beam 
constrained mass sideband data sample defined below. 

Continuum background is reduced using event shape variables. We quantify the event 
topology with modified Fox- Wolfram moments defined as 

hT = PjPkPi (cos 9jk) , 

where i enumerates B signal candidate particles {s particles) and j and k enumerate the re- 
maining particles in the event (o particles); pi is the ith particle's momentum, and Pi{cos9ij) 
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is the /th Legendre polynomial of the angle 6ij between particles i and j. The hf° terms 
contain information on the correlation between the B candidate direction and the direction 
of the rest of the event. The odd h"" terms partially reconstruct the kinematics of the other 
B in the event while the even terms quantify the sphericity of the other side of the event. 
We create a six-variable Fisher discriminant called the Super Fox- Wolfram defined as 

l=2A \"'0 J 1=1-4 \"'0 ) 

where ai and (5i are the Fisher coefficients. 

The SFW variable is combined with the B flight direction with respect to the beam 
axis, cos 6b, to form a single likelihood 

CBB = 'C{SFW)sBy<'C{cOs9B)BB 

for signal and an equivalent product for continuum, Cgg. Continuum background is sup- 
pressed by cutting on the likelihood ratio 

LR 



^BB + ^qq 



These variables are shown in Fig. |^. The signal probability density functions (PDF) are de- 
rived from Monte Carlo (MC); the continuum PDFs are taken from sideband data discussed 
below. The SFW PDFs are modeled as the sum of a simple Gaussian and an asymmetric 



Gaussian for both signal and continuum; the cosOb PDF is modeled as a second-order 
polynomial for signal and is flat for continuum. We make separate requirements on LR for 
each mode depending on the expected background determined using sideband data. As an 
example. Fig. ^ shows the B^ K'^tt' data sample before and after imposing the LR > 0.8 
requirement. 

Table | lists the reconstruction, particle identiflcation, and continuum suppression effi- 



ciencies for each final state |£8|. The reconstruction and continuum suppression efficiencies 
are determined using a GEANT-based MC [^. The error in the reconstruction efficiencies 
are determined by embedding MC generated particles into hadronic event data and compar- 
ing the efficiencies between the embedded events and the default MC and also by measuring 
the relative yields of D decays to various final states. The charged track, 7r°, and se- 
lection criteria efficiencies are tested by measuring the D event yields before and after each 
cut is applied. Further comparisons are made between kinematic distributions of particles 
in sideband data (discussed below) and continuum MC events. Based on the results of these 
studies, we assign a relative systematic error in the reconstruction efficiencies of 2.5% for 
charged tracks, 6.3% for Kg mesons and 7.3% for vr" mesons. The relative systematic error 
associated with the continuum suppression cut is 4% which is determined by taking the 
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ratio of D^ti^ yields in data after and before continuum suppression is applied and 

comparing to the MC efficiency. The -k^-k^ final state includes an additional relative 

systematic error of 10% to account for difficulties in triggering and hadronic event selection 
for this mode. 

A critical feature of the analysis is the measurement of the particle identification efficiency 
and fake rate. These are determined using nearly pure samples of K and vr mesons tagged 
using the continuum D*^ decay chain D*^ —>■ D^tt~^, —>■ K~Ti^ . Figure |^ shows the 
K~TT^ invariant mass distributions before and after applying PID cuts. For tracks in the 
B hh signal region of 2.4 GeV < p < 2.8 GeV (1.5 GeV < pi^b < 4.5 GeV) the K 
efficiency and fake rate are ex = 0.86 and fx = 0.086 (true K fakes vr); the vr efficiency 
and fake rate are e^r = 0.88 and /^r = 0.071 (true vr fakes K). The relative systematic errors 
are 2% in the efficiencies and 4% in the fake rates. These errors are directly related to the 
sample purity. Figure |^ also shows B^ —>■ K^tt^ and tt+tt" MC events before and after 
applying pion identification cuts on both tracks. 

IV. B RECONSTRUCTION AND YIELD EXTRACTION 

To reconstruct B mesons we form two quantities: the energy difference, AE = Eb — 
-Ebeam, and the beam constrained mass, m^c = \/-£'beam ~ Ps' where -Ebeam = = 5.29 

GeV, and Eb and pb are the reconstructed energy and momentum of the B candidate in 
the center of mass frame. These are shown in Fig. ^ for 5"^^+^ — > f^+vr^*^'') and tt+tt^*^'^^ 
MC. Modes containing vr'^s have a tail extending into the negative AE region due to shower 
leakage out of the back of the calorimeter and photon interactions with the material in front 
of the calorimeter. We calculate the energy of final state charged particles using a pion mass 
assumption. This shifts AE by —45 MeV for each charged kaon in the final state. The 
signal yields are extracted by a binned maximum-likelihood fit to the AE distribution in 
the region 5.271 < nibc < 5.289 GeV ( m^c > 5.270 GeV for modes containing vr^s) and —300 
MeV < AE < 500 MeV. The yields are verified by fitting in the AE signal region. A 
sideband region of 5.2 GeV < m^c < 5.26 GeV is used to study the continuum background 
in the AE distribution, while a sideband of 150 MeV < AE < 500 MeV is used to study 
the continuum background in the mf,c distribution. 

The AE fits include four components: signal, crossfeed from other mis-identified B ^ hh 
decays, continuum background, and backgrounds from multi-body and radiative charmless 
B decays. These are shown in Fig. ^. The crossfeed component is shifted from the signal 
component by 45 MeV as described above. The charmless B decay background is dominated 
by events where the B meson decays to an hhii final state such as pvr or /o(980)A' where 
one pion is not reconstructed. This shifts the charmless B background by at least the mass 



8 



of the missing pion. We expect no backgrounds from b ^ c decays based on a large MC 
sample. 

For charged particle final states, the AE signal is modeled with a Gaussian. For modes 
containing vr'^s, the signal is modeled as the sum of a primary Gaussian and a secondary 
asymmetric Gaussian. The mean positions of the two are equal and the +AE a of the 
asymmetric Gaussian is constrained to equal the a of the primary Gaussian. The crossfeed 
component has an equal shape, shifted by 45 MeV for each mis-identified particle. 

The widths of the AE signal distributions are determined using inclusive high momentum 
K~Ti^ , K'TT^n^, and Kgir'^ decays after requiring the D daughter particles 

to have a momentum range similar to B ^ hh candidate particles. These distributions are 
shown in Fig. |[ Comparisons between the D mass widths in MC and data are used to scale 
the B hh AE MC widths. This procedure is also used to determine the ratio of primary 
to secondary Gaussians for modes containing vr'^s. 

The peak positions of the AE signal Gaussians are a function of the beam energy and 
the momentum scale. The beam energy is determined using the peak position of the m^c 
distribution for the B^ D^n^, K^tt^ data sample shown in Fig. |[ The momentum 

scale is determined using the peak positions of the inclusive D mass spectra discussed above 
as well as the AE distribution for the B^ D^tt^ data sample also shown in Fig. ^. 

The continuum background is modeled with a second-order polynomial with coefficients 
determined from sideband data. Figure demonstrates the validity of this method by 
comparing the continuum AE background shape in on-resonance mf,c sideband data to the 
shape in the rribc signal region in the off-resonance data sample. Backgrounds from charmless 
B decays are modeled by a smoothed MC histogram. 

For all final states except K'^n^ and 7r"*"7r°, the normalizations of the four components are 
the only free parameters in the fits. The significance of the signal yield above background 
is determined by re-fitting the AE distribution without a signal component and comparing 
the maximum likelihoods of the two fits. Due to the large overlap of the signal and crossfeed 
components in the K~^7t^ and vr^^vr'^ signals, we perform a simultaneous fit to the K'^n'^ and 
tt+tt" AE distributions constraining the crossfeed to the expected values based on the PID 
fake rates. 

The rribc distribution provides no discrimination among the three B decay components. 
The sum of the three components is modeled with the same functional form as the AE 
signal shapes discussed above. We parameterize the continuum background with a function 
that behaves like phase space near the endpoint (the ARGUS shape [p^ ). 
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V. RESULTS 



Figures and TU show the AE and m^c distributions for the Kn, tttt and KK 

final states, respectively. The AE signal yields and the significance above background are 
listed in Table |I| [0. Using these results and the efficiencies listed in Table |I|, we derive 
the branching fractions listed in Table |l| based on the data sample of 31.7 million BB 
events In all cases, the m^c fits give consistent results. The systematic error in the 

fitting procedure is determined by varying the parameters of the fitting functions within 
their errors and measuring the change in the signal yield. The deviations from the nominal 
yields are typically 1 to 2 events. These deviations, along with the error in the efficiencies 
and N{BB) (1%) are added in quadrature to give the systematic error in the branching 
fractions. For modes with significance below Sex B^, we report 90% confidence level upper 



limits calculated with the efficiency and N{BB) reduced by their systematic errors. 

For the K+n^/n+n^ simultaneous fit, we re-fit the distributions after removing the con- 
straints. The central values of the signal yields differ by at most 2.3 events. These deviations 
are also included in the systematic error. 

The effects of backgrounds from charmless B decayscan be demonstrated by re-fitting 
the AE distributions in the region AE > —130 MeV without a charmless B background 
component. For modes with vr^s in the final state, the yields deviate by as much as 12%, 
clearly indicating the need for these components in the fit. There is almost no deviation 
in the yields of final states that do not include a 7r°. In each fit, the measured charmless 
B background yield agrees with the expected values based on independent measurements 
of these modes. However, the errors on the fitted yields and the expected values are both 
large. 



In Table III, we list ratios of partial widths among the Kn and tttt modes 31 . The 



correlations between the numerator and denominator are included in the systematic error 
calculation. The systematic error includes a 2.5% fractional error from the ratio of charged 
to neutral B meson lifetimes where we have used /t^ = 1.091 ± 0.027 to convert the 
ratio of branching fractions to the ratio of partial widths. 

We measure partial-rate asymmetries in all measured fiavor-specific modes. The AE 
distributions are shown separately for B and B modes for the K'^n^, K^ir'^, and n^ir^ 
final states in Fig. |Tl|. Figure [l^ shows the corresponding distributions for the K^n^ final 
states. The fitting results and partial-rate asymmetries are listed in Table Here, the 90% 
confidence intervals assume Gaussian statistics and are expanded linearly by the systematic 
error. The systematic errors are dominated by fitting systematics but also include a 1% 
contribution, added in quadrature, to account for possible detector-based asymmetries, as 
discussed below. In the K'^n^ final states, the asymmetry is corrected by a dilution factor 
of 0.984 ± 0.001, due to double misidentification of K~^tt~ as K~7t~^. 
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Four samples are used to verify the symmetric performance of the Belle detector for high 
momentum particles. An inclusive sample of tracks in the two-body decay momentum bin 
2.4 GeV< p < 2.85 GeV is used for tracking efficiency tests before and after PID cuts are 
applied. Events in m^c sideband data further test the reconstruction efficiency along with 
the continuum suppression cut efficiency. We also check the difference between inclusive 
high momentum D^^~^^ K~tt^, Kgir^, and K^n^n^ decays and their charge conjugates 
that test for asymmetries in the detector resolution. The entire reconstruction procedure is 
applied to the B Dvr^; D K^n^, K^tt^, and iT^ n'^ir'^ data samples. The results are 
listed in Table |V]. The inclusive track sample yields an asymmetry in the track reconstruction 
efficiency of {N{h-) - N{h+))/ {N{h-) + N{h+)) = (-3.6 ± 0.3) x IQ-^. Considering the 
statistical precision of the current data set, we ignore this very small asymmetry. We also see 
a 1.8 MeV shift between the inclusive D~ — >■ K^7r~ and — > Kg7r~^ mass peaks indicating 
an approximate 0.1% momentum scale difference between positive and negative tracks. This 
shift has been taken into account when determining the error in the peak position of the 
AE distributions for modes with odd numbers of tracks in the final state. Furthermore, 
as shown in the Table 0, the shift does not result in an asymmetry between the efficiency 
to reconstruct the two flavors. All other mass and width parameters are consistent within 



the errors between the two flavors for the three D decay channels. Figure [T^ shows the 
samples most relevant to the —* K^tt^ flnal states. With the exception of the inclusive 
track sample, all are consistent with zero asymmetry within the error and we conclude that 
detector based asymmetries are below 1%. 



VI. DISCUSSION AND CONCLUSIONS 

We have presented measurements of the branching fractions of all B — » Kit flnal states 
and the B tt^tt^ and tt^tt*^ flnal states. We see no signiflcant evidence for the decays 
B^ tt^tt" ot B —>■ KK and set 90% confldence level upper limits on their branching 
fractions. Furthermore, we see no clear evidence for partial-rate asymmetries between the 
B and B decay amplitudes in these modes. 

The partial-rate asymmetry between B^ K^7i~ and 5^ K^-k^ of 0.46 ±0.15 ±0.02 
has a non-zero signiflcance of 2.9cr. Since this is below 3(j, we defer claiming evidence for an 
asymmetry until a larger data sample is collected. 

The r(7r^7r^)/r(i^'^7r^) and r(7r^7r^)/2r(7r^7r°) partial width ratios are signiflcantly 
below 1. This could be an indication of large destructive interference between Tree and 
Penguin amplitudes for the B^ 7i^7t~ decay although theoretical uncertainties pertaining 
to the relative sizes of the interfering amplitudes are still large 

Knowledge of the branching fraction for 5° — > vr'^vr^ is required for an isospin analysis of 
B ^ TXTX decays and for the extraction of the CP violation parameter sin 202 |]Tn[- With the 
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present data set we set a limit of B{B^ — > 7r°7r°) < 6.4 x 10 ^ at the 90% confidence level. If 
the excess of ~ 13 events is indeed due to S° — > 7r°7r° events, the corresponding branching 
fraction would be B(B^ n^n^) = (3.2 ± 1.5 ± 0.7) x 10"^ 
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FIG. 1. The TT+TT"" (left) and 77 (right) mass spectra for 5+ — Kgh'^ and /i+7r° candidates 
in the beam constrained mass sideband data sample. The 7r+7r~ distribution is modeled as the 
sum of two Gaussians for true Kg candidates while the background is modeled as a second-order 
polynomial. The weighted average resolution of the two Gaussians is 3.4 MeV for the Kg mass 
peak. The 77 distribution is modeled as the sum of a primary symmetric Gaussian and a secondary 
asymmetric Gaussian [24] for vr*^ candidates and a second-order polynomial for background. The 
weighted average resolution of the two Gaussians is 8.9 MeV for the ir^ mass peak. In both 
distributions the solid curves are the sum of signal and background components while the dashed 
curve is the background component. 
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FIG. 2. Continuum suppression variables: SFW (top left), cos 6b (top right), and the combined 
likelihood ratio (bottom left). The solid curves are the signal PDFs derived from MC. The dashed 
curves are the continuum PDFs derived from sideband data. The open points are the S+ D°7r+, 
D° K~7r+ data sample. The solid points are off-resonance data. The bottom right distribution 
is the beam constrained mass distribution for the K'^it~ data sample before and after 

requiring the likelihood ratio cut LR > 0.8. 
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FIG. 3. The D*+ D°7r+, D° K~7r+ decay sample used to determine the PID efficiencies 
and fake rates, shown in the top row. The distribution on the left is the D mass distribution where 
the kaon daughter is required to have a similar momentum/cos^ distribution as the two-body 
B ^ hh decay daughters. The open histogram contains all candidates. The hatched histogram 
contains events where the kaon passes PID cuts. The solid histogram contains events where the 
kaon is mis-identified as a pion. The distribution on the right is the corresponding figure for the 
pion daughter. The bottom row shows the combined —> K~^tt~ and 7r"*"7r~ MC AE distribution 
assuming a 4:1 K~^TT~:iT~^7r~ production ratio. The distribution on the left is before PID cuts are 
applied. In the distribution on the right, both tracks are required to be identified as pions. The 
solid histogram is the sum of K^n^ and tt+tt", the dashed histogram is the K~^Tr~ component, 
and the hatched histogram is the ■K^7r~ component. 
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FIG. 4. Kinematic reconstruction variables (left) and AE (right) for K~^Tr~ 

and 7r"'"7r~ MC (top) and K^n^ and tt+tt^ MC (bottom). The solid histograms are 

TT+TT ^ ' MC. The dotted histograms are ofF-resonance data. The dashed histograms 
are for B^^'^^ K'^Tr~^^^ MC events, which are indistinguishable in the m^c distribution but 
shifted by —45 MeV in AE due to pion mass assignment to the kaon track. The dot-dashed 
histograms represent background from multi-body charmless B meson decays. 
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FIG. 5. Mass spectra for inclusive K~7r^ (top left), — > Kgir^ (top right), and 

D° — s- K~'k'^it^ (bottom) used to determine AE fit parameters. In each case, momentum cuts are 
placed on the D daughter particles to simulate the momentum of B daughter particles. 
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FIG. 6. The S+ — D°7r+, K~7r+ data sample used to determine rribc (left) and AE 

(right) fit parameters. The AE distribution contains backgrounds from B° D°7r+7r~. 
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FIG. 7. The AE distribution for 5+ — Kg7r~^ candidates in rubc sideband data 
(dashed-histogram) compared to the same distribution in the nibc signal window for off-resonance 
data (black circles). The solid curve is a second-order polynomial used to parameterize the shape. 
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TABLE I. 


Efficiencies to reconstruct the B 


hh modes. 


Listed are the efficiencies for re- 


construction, 


particle identification, 


continuum suppression, the final combined efficiencies, and 


products of efficiency times intermediate branching 


fraction. 






Mode 




Rec. 


PID 


qq 


Final 


ex B.F. 






0.73 


0.76 


0.55 


0.31 


0.31 






0.75 


0.77 


0.52 


0.30 


0.30 


K+K- 




0.71 


0.74 


0.38 


0.20 


0.20 


i\ vr 




n /I Q 
U.4o 


U.OD 


0.39 


0.14 


0.14 






0.46 


0.88 


u.oy 


U. ID 


U. ID 


K^TT+ 




0.53 


0.88 


0.68 


0.32 


0.11 


K'>K+ 




0.51 


0.86 


0.38 


0.17 


0.06 






0.34 


1 


0.69 


0.23 


0.08 






0.31 


1 


0.41 


0.13 


0.13 




28| 


0.37 


1 


0.54 


0.20 


0.04 
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FIG. 8. The /S.E (left) and (right) fits to the B Ktt event samples. The sum of the 
signal and background functions is shown as a solid curve. For the AE distributions, the dashed 
curve represents the signal component, the dotted curve represents the continuum background, 
and the hatched histogram represents the charmless B background component. For the K^tc~ 
and K~^TT^ distributions, the crossfeed components from ir^ir^ and Tr^Tr" are shown by dot-dashed 
curves centered 45 MeV above the signal components. For the nn,c distributions, the continuum 
background is represented by the dotted curve while the sum of signal, charmless B background, 
and crossfeed components is shown by the dashed curve. 
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FIG. 9. The AE (left) and (right) fits to the B tttt event samples. The sum of the 
signal and background functions is shown as a solid curve. For the AE distributions, the dashed 
curve represents the signal component, the dotted curve represents the continuum background, 
and the hatched histogram represents the charmless B background component. For the tt+tt" and 
7r+7r° distributions, the crossfeed components from K'^ir~ and K'^tt'^ are shown by dot-dashed 
curves centered 45 MeV below the signal components. For the m^c distributions, the continuum 
background is represented by the dotted curve while the sum of signal, charmless B background, 
and crossfeed components is shown by the dashed curve. 
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FIG. 10. The AE fits to the B KK event samples. The sum of the signal and background 
functions is shown as a sohd curve, the dotted curve represents the continuum background. In 
the K^K~ distribution, the dot dashed curve represents the K+tt" crossfeed. In the 

5+ K'gK'^ distribution, the hatched histogram represents the charmless B background. In the 
B^ KgKg distribution, the dashed curve represents the signal component. 



TABLE II. 


Signal yields, significance 


above background, and branching 


fractions for B ^ hh 


modes assuming 


I equal production fractions for neutral and charged B meson pairs. We report 90% 


confidence level 


upper limits for the B^ — > tt'^tt'^ and B KK decays. 






Yield 


Sig. 


B.F. (xlO-5) 


K+TT- 


217.6ti?J 


16.4 


2.25 ± 0.19 ± 0.18 




OO-O_10.7 


6.4 


1.30ig:^^ ± 0-13 


K\+ 


D'-'-'-lO.l 


7.6 


1.94to|i±0.16 






2.8 


0.80lg|? ± 0.16 




n+11.6 
^^■^-10.9 


5.4 


0.54 ±0.12 ±0.05 


tt+ttO 


36.7ll^J 


3.5 


0.74l°:i ± 0.09 






2.4 


< 0.64 


K+K- 







< 0.09 


K^K+ 







< 0.20 




9+2-9 

'-'•^-0.9 





< 0.41 
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FIG. 11. The AE distributions for B (left) and B (right) candidate decays to flavor-specific 
final states. The sum of the signal and background functions is shown as a solid curve. The dashed 
curve represents the signal component, the dotted curve represents the continuum background, 
and the hatched histogram represents the charmless B background component. The crossfeed 
components are shown by dot-dashed curves centered 45 MeV from the signal components. 
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FIG. 12. Distributions for ^s^r^. The top row contains AE distribution for B (left) 

and B (right) candidates. The corresponding rribc distributions are shown in the bottom row. The 
sum of the signal and background functions is shown as a solid curve. The dashed curve represents 
the signal component, the dotted curve represents the continuum background, and the hatched 
histogram represents the charmless B background component. 



TABLE III. Ratios of partial widths among the various B Kit and vrvr final states assuming 
equal production fractions for neutral and charged B meson pairs. The ratios of branching fractions 
are converted to ratios of partial widths using t~^/t^ = 1.091 =b 0.027 [34]. 

r2/ri 

TT+TT-fK+TT- 0.24+0^1 ± 0.02 

2K+7r^/K^n+ 1.34 ± 0.33 toil 

{T+/r^)K+7T-/K''7r+ l-27l°-| ± 0.10 

(r+/r°)7r+7r-/27r+7r0 0.40 ± 0.15 ± 0.05 

(r+/T0)7r07r77r+7r0 < 0.83 (90% C.L.) 



26 



300 



J" 

!§ 200 



100 



-0,25 





1 1 1 1 1 1 1 1 1 








0.00 0,25 

AE (GeV) 



300 



0,50 




0,00 0,25 

AE (GeV) 



0,50 



800- 
600- 
400- 




1,78 1,82 1.86 1,90 
m(K^%~)(GeV) 





800- 


'- ■> 

- 1 




600- 










\ 


400 1 


Evem 


200- 











V 







1.78 1.82 



1.86 1.90 



1.94 



m(K^%)(GeV) 

FIG. 13. Control samples for the ^s^^ -Acp measurement. The top row contains the 

distributions for B" D^-k' , K'^tt" on the left and the charge conjugate decay on 

the right after applying identical PID and continuum suppression cuts as for the K^tt^ analysis. 
The bottom row contains inclusive D~ — K^tt" (left) and — K^tt'^ (right) mass spectra after 
making momentum cuts on the D daughter particles to simulate the momentum of B daughter 
particles. 

TABLE IV. Partial-rate asymmetries. Listed are the number of signal events for each final 
state, the .Acp values with errors, and their 90% confidence intervals, listed on the following line. In 
the K^TT^ final states, the asymmetry is corrected for the dilution due to double mis-identification. 



Mode 



NiB) 



NiB) 



Acp (90% C.L.) 



102.8 ± 12.6 
28.7 ±7.8 
49.5 ± 8.4 
24.2 ± 8.4 



115.0 ±13.3 
30.1 ± 7.7 
18.6 ±6.3 
13.0 ± 7.3 



-0.06 ± 0.09l|j:|j2 

-0.21 : 0.09 
-0.02 ±0.19 ±0.02 

-0.35 : 0.30 
0.46 ± 0.15 ± 0.02 
0.19 : 0.72 
0.30 ± 0.30lg:g^ 
-0.23 : 0.86 
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TABLE V. Tests of detector based asymmetries. The first row is based on an inclusive track 
sample with 2.4 GeV < p < 2.85 GeV in the center of mass frame. The following two rows are the 
asymmetries in the same sample after applying particle ID. The fourth row is the asymmetry for 
inclusive D meson decays to high momentum K+tt", K+7r~7r°, and KgTr~^ final states. The fifth 
and sixth rows are the asymmetries in the m^c sideband and the —>■ D^^~^Tr~^ data sample 

before and after the LR cut is applied. 



test sample 


Acp 




high p tracks 


(-3.6 ± 0.3) X 10-3 




with K pid 


(-3.2 ±0.5) X 10-3 




with TT pid 


(-3.7 ±0.3) X 10-3 




high p D decays 


(-2 ± 3) X 10-3 






before LR cut 


after 


rrihc sideband 


(-0.07 ±0.17) X 10-3 


0.01 ± 0.05 


B — Dn^ 


-0.045 ± 0.025 


-0.055 ± 0.027 
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